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ABSTRACT 


We discuss the observational evidences of the morphological transformation of Spirals into 
SO galaxies in the cluster environment exploiting two big databases of galaxy clusters: WINGS 
(0.04 < < 0.07) and EDisCS (0.4 < 2 ; < 0.8). 

The most important results are: 1) the average number of SO galaxies in clusters is almost a 
factor of ~ 3 - 4 larger today than at redshift z ~ 1; 2) the fraction of SO’s to Spirals increases 
on average by a factor ~ 2 every Gyr; 3) the average rate of transformation for Spirals (not 
considering the infall of new galaxies from the cosmic web) is: ~ 5 Sp^SO’s per Gyr and ~ 2 
Sp-^E’s per Gyr; 4) there are evidences that the interstellar gas of Spirals is stripped by an hot 
intergalactic medium; 5) there are also indirect hints that major/minor merging events have played 
a role in the transformation of Spiral galaxies. 

In particular, we show that: 1) the ratio between the number of SO’s and Spirals {Nso/Nsp) in 
the WINGS clusters is correlated with their X-ray luminosity Lx', 2) that the brightest and massive 
SO’s are always close to the cluster center; 3) that the mean Sersic index of SO’s is always larger 
than that of Spirals (and lower than E’s) for galaxy stellar masses above 4) that the 

number of E’s in clusters cannot be constant; 5) that the largest difference between the mean 
mass of SO’s and E’s with respect to Spirals is observed in clusters with low velocity dispersion. 

Finally, by comparing the properties of the various morphological types for galaxies in clusters 
and in the field, we find that the most significant effect of the environment is the stripping of the 
outer galaxy regions, resulting in a systematic difference in effective radius and Sersic index. 

Keywords: clusters of galaxies; galaxy formation; galaxy morphology 


1 INTRODUCTION 


The morphological class of SO galaxies has been largely debated since [Hubble ( |1936[ ) supposed their 
existence in “Realm of the Nebulae”, and Sandage ( 1961 ) placed them as transition types in the tuning 
fork diagram. [Spitzer & Baade | ( [r951| ) first suggested that SO’s might originate from Spirals stripped out of 
gas and that this class of objects could form a dustless sequence that parallels the Sa-Sb-Sc sequence in the 
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Hubble tuning fork. This idea was further developed by van den Bergh (1976) who proposed the class of 
“anemic Spirals” as intermediate objects between normal SO’s and Spirals. 


Historically there are essentially two main views concerning the origin of SO galaxies: 1) SO’s are the end 
product of the transformation of Spiral galaxies removed of their gas content by means of gravitational or 
non gravitational forces, or by secular evolution; 2) the SO class is intrinsic, i.e. formed by objects that 
differ from Spirals since their formation. 


The proponents of the first view found support in the following observational facts: 1) a significant 
fraction of today’s early-type galaxies likely evolved from later-type galaxies at relatively recent cosmic 
epochs. Hubble Space Telescope (HST) observations revealed that Spirals are proportionally much more 
common (a factor of ~ 2 — 3), and SO galaxies much rarer, in distant than in nearby clusters (see e.g. 


Dressier et al. 1997; Fasano et al. 2000t [Treu et al. [ [2003t [Postman et al. [ |2005t [Smith et al. [ 2005; 

Desai et al. 

2007); 2) several evidences exist that Spirals are stripped of their HI gas in dense environments 

(see e.g. Davies & Lewis 1973; Haynes et al. [[1984;[Giovanelli & Haynes [ 1985 

Solanes et al. 2001; 

Koopmann & Kenney [ 2004 

Growl et al. [ 2005 1 Chung et al. 2007[ Sun et al. 

20071 Yagi et al. 2007 

Rasmussen et al. [ 2008t [Vollemer et al. [ |2009t [Yoshida et al. 2008 Smith et al. [|2010t[Sivanandam et al 


[ 2010). Loosing their gas Spiral galaxies rapidly transform into SO’s owing to the quenching of their star 
formation (SF) (see e.g. [Boselli & Gavazzi \ |2006t [Book & Benson |2010[ |Fabello et al. [ |2011[ [Goncalves 
et al. [ 20 12[ [Mendel et al7 2013[[Bekki 2014[ Hirschmann et al. 2014 Wu et al. 2014); 3) the specific 


frequency of globular clusters (GCs) in present day SO’s is consistent with the idea that these galaxies are 


formed when the gas in normal Spirals is removed (Barra et al. 2007); 4) the observed mass-metallicity 


trends in SO galaxies are consistent with a merging formation of these galaxies (Chamberlain et al.^ 2012); 


6) the stellar kinematics of SO galaxies inferred from Planetary Nebulae (PN) appears consistent with 
those of Spirals, indicating that they could evolve directly from such systems via gas stripping or secular 


evolution (Cortesi et al. I 2011). 


On the other side the proponents of the second view claim that: 1) the bulge-to-disk ratios (B/D) of the 
two morphological classes are too different for being explained by the simple occurrence of gas stripping 
(see e.g. [Burstein [|1979} . If SO’s are evolved Spirals, their bulge luminosities must have been physically 
enhanced ( Christlein & Zabludoff [[2004 ); 2) many edge-on SO’s are characterized by thick disks that 
are absent in the control samples of Spirals; 3) the relation between density and morphology is a weak 
function of density, so that a significant percentage of SO’s are in regions where the density and temperature 


is too low for a significant stripping of the gas (Dressier 1980); 4) the large ages and magnesium-to 


iron overabundance of the outer disks of nearby SO’s contradict the commonly accepted scenario of 
SO formation from Spirals by quenching their star formation during infall into dense environments at 
intermediate redshifts, z = 0.4 — 0.5 (4-5 Gyr ago) ( [Silchenko et al. [[2012[ ); 5) the evolution in the 
number density of SO’s from 2 ; = 0.4 — 0.5 to 2 ; = 0 is the result of morphological classification errors 
( Andreorfl 1998); 6) neither the median ellipticity, nor the shape of the ellipticity distribution of cluster 
early-type galaxies evolve with redshift from 2; ~ 0 to 2; > 1, i.e. over the last ~ 8 Gyr ([Holden et al. [ 


2009[), suggesting that the bulge-to-disk ratio distribution for cluster early-type galaxies does not change 


with time. 


The controversy on the SO’ origin is not yet settled after more than 70 years! There is however today a 
wide consensus that more than one process can operate the transformation of Spirals into SO’s and that a 
significant role is played by the environment. 


This is a provisional file, not the final typeset article 
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In low density environments the transformation of Spirals into SO’s eould be the end produet of 
seeular evolution. This idea is supported by the faet that many SO’s in this environment host pseudo¬ 
bulges (Laurikainen et al. 2006; Graham [ 2013), i.e. bulges with bulge-to-total luminosity ratios and 
eoneentrations typieal of late-type spiral galaxies (B/T < 0.2 and n ~ 1). These bulges often eontain 
embedded disks, inner spiral patterns, nuelear bars, star formation, and present rotational support similar 


to that of Spiral bulges (see Kormendy & Kennieutt 2004). These properties are usually attributed to 


the evolution indueed by bars, whieh send material towards the eenter of the galaxy, enhaneing the bulge 


eomponent (Pfenniger & Norman [ 19901, even if often these bulges do not seem to have the eoneentrations 


and sizes of typieal SO’s bulges (see g.g. |Eliehe-Moral et al. [|2012[|2013[). However,[Pfenniger & Friedll 


[ (1993) found that the buekling instability ean lead to vertieally thiek inner bar eomponents, whieh are 
several faetors larger than the nuelear star forming struetures. They have similar sizes as the boxy/peanut 
bulges in the SOs (and in other Hubble types) in the edge-on view. Furthermore, seeular evolution offers a 
viable meehanism to produee the thiek disk eomponent observed in many edge-on SO’s. 

In the eluster environment several meehanisms ean operate the quenehing of SF and start the transition of 
Spirals toward the SO elass: 1) gas loss triggered by: galaetie winds as a eonsequenee of star formation 
or AGN aetivity (see, e.g. [Veilleux et al. [ |2005[ |Ho et al. 1 12014[ [Fogarty et al. |, |2012|), ram pressure 


stripping, i.e. the interaetion between the interstellar medium (ISM) of the galaxy and the intergalaetie 


medium (IGM) (Gunn & Gott 1972), strangulation, i.e. the removal of the hot gas halo surrounding the 


galaxy via ram pressure or tidal stripping by the halo potential ( jFarson et al. [|1980tlBalogh et al. [|2000| ), 


thermal evaporation (Cowie & Songaila , 1977), and turbulent/viseous stripping (Nulsen [ 1982). 2) strong 


tidal interaetions and mergers, although more frequent in groups than in elusters ( [Barnes & Hernquist [ 
1992; Mazzei et al. 2014a[[b ), tidal effeets from the eluster as a whole ( Byrd & Valtonen1[1990| ), and 


harassment, i.e. the eumulative effeet of several weak and fast tidal eneounters (Moore et al. , 1996). 


Ram pressure is one of the most aeeredited aetors of the transformation of Spirals in elusters. Its 
importanee resides in the faet that only the gaseous eomponent of disk galaxies is affeeted. The ram 
pressure exerted by the IGM is a key ingredient for understanding the star formation history (SFH) of 
galaxies in elusters (see e.g., Bekki [ 2014), beeause temporal and spatial variations of the SF in disks 
ean be produeed by ram pressure (see e.g., Koopmann & Kenney [[2004). Observations suggest that the 


Ha-to-optieal disk-size ratio is smaller for Spirals with higher degrees of Hl-defieieney (Fossati et al. 
2m3| ). 


The merging hypothesis is also quite popular and old (see e.g., Biermann & Tinsley ( [1975 ); meehanism 
defined at that time ’’eollision between galaxies”) and often advoeated for explaining the existenee of SO 
galaxies in low density environments (e.g. groups and held) where the meehanisms typieal of elusters are 
less effeetive. This possibility found support in the diseovery of ’’Polar Ring SO’s”, i.e. lentieular galaxies 
hosting struetures likely due to a seeond event sueh as the transfer of mass from a eompanion galaxy 
during a elose eneounter as well as a true merger (Sehweizer et al. [ 1983). Bekki [([1997 ) stressed that a 


partieular orbit eonhguration eould transform two late-type spirals into one SO galaxy with polar rings (see 
also Mapelli et al. [ 2015[ ). A related phenomenology, namely SO’s hosting an unexpeeted gas eomponent 
eounter-rotating or strongly kinematieally deeoupled with respeet to the stellar body of the galaxy, has also 
been deteeted ( Bertola et al. [ 1992). Even more, numerieal simulations by [Bekki (1998) indieate that the 
major meehanism for the SO ereation is galaxy merging between two spirals with unequal mass; if this is 
true, unequal-mass galaxy mergers provide the evolutionary link between the large number of blue spirals 
observed in intermediate-redshift elusters and the red SO’s in the present-day ones. 
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An additional piece of evidence favoring the merging idea comes from the correlation observed in 
very nearby galaxies between UBV color residuals and the Schweizer’s hne-structure S index, especially 
because the detected systematic variations are not limited to the nuclei, but occur globally in the stellar 


populations of galaxies (Schweizer & Seitzer 1992). Again with reference to observed colors, the so-called 


blue E/SO galaxies studied by |Aguerri et al. | ( |2010p , i.e. objects having a clear early-type morphology on 
HST/ACS images but with a blue rest-frame color, do resemble merger remnants probably migrating to 
the red-sequence in a short time-scale, when their mass is comparable (or higher) than that of the Milky 


Way. Analogously Prieto et al. (2013) identify several pieces of evidence favoring the same hypothesis 
i.e. that major mergers have played a dominant role in the definitive buildup of present-day E-SO’s with 

M* > IO^^Mq at 0.6 < z < 1.2. 

At any rate the frequency of other photometric hints indicating a major role of merging events in SO 
galaxies formation is still controversial. Eor instance the claim by |Borlaff et al. (2014l, who find that—in 
the context of their N-body simulations—nearly 70% of their SO-like host anti-truncated (Type III) stellar 


disks does not find support in the observational works by Erwin et al. (2008) as well as by Gutierrez 


et al. (2011|), who at variance conclude that the difference in the fractions of Type-III profiles between 


early (presumably SO-SOa) and late-type spirals is not statistically significant, in other words, it should not 


reasonable to conclude that Type-III photometric profiles are a special characteristic of SO’s. Maltby et 


al. (2015) reports for a sample of ~ 280 field and cluster SO’s similar fractions of Type profiles in both 


environments. 


On the other side Querejeta et al. (2015) have shown that major mergers among binary galaxies can 
rebuild a bulge-disk coupling in the SO remnants in less than ~3 Gyr, after having destroyed the structures 
of the progenitors, whereas minor mergers directly preserve them. The resulting structural parameters span 
a wide range of values, and their distribution is consistent with observations. In their simulations many 
remnants have bulge Sersic indices ranging 1 < n < 2, flat appearance, and contain residual star formation 


in embedded discs, a result which agrees with the presence of pseudo-bulges in real SO’s. Querejeta et 


al. (2015b I also find that major mergers can explain the difference in the stellar angular momentum and 


concentration between Spirals and SO’s. 


Disentangling the mechanisms at work in the transformation of Spirals into SO galaxies is not easy 
even in very nearby objects. Eor more distant galaxies a combination of efforts are required, both in the 
sense of increasing the observational data available for galaxies in clusters at different redshifts, and in 
producing numerical simulations of increasing complexity. With this paper we want to contribute to this 
effort, discussing the data coming from two clusters surveys, WINGS (Wide field Nearby Galaxy cluster 
Survey ; |Fasano et aLl|2006| ) at redshifts 0.04 < 2 ; < 0.07 and EDisCS (the ESO Distant Cluster Survey; 
White et aT] |2005| ) at 0.4 < z < 0.8, that are particularly important for their ample databases. The idea 


behind this work is that by comparing the mean properties of galaxies of different morphological classes 
per bin of masses, cluster-centric distance, and redshift, we can find a statistically significant proof for the 
origin of the SO class. 


We start with Sec. and its subsections discussing the WINGS and EDisCS data samples used for the 
present analysis, as well as the characteristics of three additional subsamples that have been used to support 
our conclusions. In Sec. we provide the main observational evidences in favor of the hypothesis that 
Spirals in clusters transform into SO’s. In Sec.|^we present a detailed analysis of the stacked WINGS and 
EDisCS data samples discussing their implication for the origin of the SO class. In Sec. we discuss the 
evidences of ram pressure stripping and evaporation, and in Sec. |^those of merging. In Sec.j^we compare 
the properties of field and cluster SO’s and, finally in Sec. we present our conclusions. 


This is a provisional file, not the final typeset article 
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Throughout the paper we have adopted the following cosmology: Hq = 70 kms ^Mpc \ VLm = 0.3 
and = 0.7. 


2 THE DATA SAMPLES OF GALAXIES IN CLUSTERS 


2.1 WINGS 


WINGS is a long-term project dedicated to the characterization of the photometric and spectroscopic 


properties of galaxies in nearby clusters. The core of the survey is WINGS-OPT (Varela et al. 2009) 


which is a set of B and V band images of a complete, X-ray selected sample of 77 clusters with redshift 
0.04 < z < 0.07. The images have been taken with the Wide Field Camera (WFC, 34' x 34') at the INT-2.5 
m telescope in La Palma (Canary Islands, Spain) and with the Wide Field Imager (WFI, 34' x 33') at the 
MPG/ESO 2.2 m telescope in La Silla (Chile). The optical photometric catalogs have been obtained using 
SExtractor and are 90% complete at V~2L7, which translates into My -|- 6 at the mean redshift of 
the survey ([Varela et al. [|2009]). The WINGS-OPT catalogs contain ~400,000 galaxies in both the V- and 


B-band. According to Varela et al. (2009 Table D.2), in the whole cluster sample the surface brightness 
limits at ISa^kg (o'bkg is the standard deviation per pixel of the background) span the ranges 24.7 -E 26.1 
(average value: 25.71) in the V band and 25.4 -E 26.9 (average value: 26.39) in the B band. SExtractor 


catalogs have also been obtained for the near-infrared follow-up of the survey (WINGS-NIR; Valentinuzzi 


et al. 


20091, which consists of J- and K-band imaging of a subsample of 28 clusters of the WINGS-OPT 
sample, taken with the WFCAM camera at the UKIRT telescope. Each mosaic image covers ~0.79deg . 
With the SExtractor analysis the 90% detection rate limit for galaxies is reached at 1=20.5 and K=19.4. 
The WINGS-NIR catalogs contain ~490,000 in the K band and ~260,000 galaxies in the J band. The 
photometric depth of the WINGS-NIR imaging is slightly worse than that of the WINGS-OPT imaging. 
Thus, for the K-band the surface brightness limit at l.Satikg spans the range 20.6 -= 21.5 (see Table 4 in 
Valentinuzzi et al. ] 2009 [ average value: 21.15). The U-band photometry has been realized up to now for 
a subsample (17 clusters) of galaxies of the WINGS-OPT sample (WINGS-UV; jOmizzolo et aI7]|2014| ), 
but new u band imaging of the WINGS clusters have been recently acquired at the VST telescope within 
OMEGAWINGS project, that is an extention of WINGS. The area covered is ~ 1 square degree and 
allows to reach up to 2.5 cluster virial radii. OMEGAWINGS is based on two OmegaCAM@VST GTO 
programs for 46 WINGS clusters. The B and V campaigns of observations are already completed, while 
the M-band programme is still ongoing. The B and V data of OMEGAWINGS are presented in Gullieuszik 
et al. (2015). 

In addition to the aperture photometry catalogs, the WINGS database includes the surface photometry 


catalogs realized with GASPHOT (D’Onofrio et al. [|2014| ). This automatic software provided for each 
galaxy and each band (with an area larger than 200 pixels) the total luminosity, the effective radius Re, the 
effective surface brightness (/i)e, the Sersic index n, and the minor to major axis ratio b/a. In addition the 
automatic tool MORPHOT (jPasano et al. 2012) produced the accurate morphological classification T of 


^40,000 galaxies. Eor more information about the WINGS database see Moretti et al. (2014). 


The medium-resolution multifiber spectroscopy of 48 clusters is another valuable property of the WINGS 


database (WINGS-SPE; Cava et al. 2009 >. Redshifts and membership were measured (see for details. Cava 


et al. [|200^ , as well as star formation histories, stellar masses, ages (see for details, [Fritz et al. [[2007[|201 f| ), 
equivalent widths, and line-indices ( Fritz et aI7] 2014[ ) for ~6,000 galaxies. Within OMEGAWINGS new 
spectra have been obtained with AAOmega@AAT using the OmegaCAM fields. So far, we have secured 
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high quality spectra for ~ 30 OMEGAWINGS clusters, reaching very high spectroscopic completeness 
levels for galaxies brighter than V = 20 from the cluster cores to their periphery (Moretti et al. in prep.). 

2.2 The WINGS stacked sample 

Using all the galaxies members of the 48 WINGS clusters with available spectroscopjQwe have built 
a stacked sample aiming at increasing the statistical significance of the results obtained for the different 
morphological classes. 

We extracted from the WINGS database the following data: the identification name (ID), the cluster name, 
the redshift (z), the scale of the images (kpc/arcsec), the value of Mpc obtained for the cluster, the 

morphological type (T), the luminosity distance (i^/um)^ the equatorial coordinates of the galaxies (RA 
and DEC), the B and V GASPHOT total magnitudes, the circularized effective radius (i?g) in kpc unit, the 
mean effective surface brightness ((/r)e) in mag arcsec“^, the distance from the brightest cluster member 
(Dbcg) in R200 unit, the star formation rate in the last 2 x 10^ Myr (SF R[MQyr~^]), the axis ratio in the 
V band (b/a), the total stellar mass (M*[0]), the luminosity weighted age {wAge[yr]), the local density 
(LD), the galactic extinction in the B and V bands, the Sersic index (n) in the V band, and the average 
B—V color. 

In total we got 2869 galaxies members of the clusters from a total number of objects of the spectroscopic 
catalog of ~ 6000 galaxies. According to their morphological classification we have 803 Elliptical galaxies 
(T < —5), 1347 SO galaxies (—4.5 < T < 0), and 719 (T > 0) Spiral galaxies. They span an interval of 
masses from ~ 10® to ~ 10^^ Mq, and an interval of distances from the center from 0 to ~ 1 i? 200 - 

The stacked WINGS sample is aimed at creating a “pseudo-cluster” having the mean properties of 
our clusters, but with many more objects at each distance Dbcg/R 200 from the center. This is possible 
because the WINGS clusters share well defined properties in terms e.g. of X-ray luminosities, velocity 
dispersions, dimensions, richness, ages, etc. This homogeneity of properties is also accompanied by a set 
of homogeneous observations and data reductions, so that the final sample is very robust from a statistical 
point of view. We can say that this pseudo-cluster provides us the mean properties of the WINGS clusters 
and its variance. This stacked sample permits the analysis of the mean properties of the morphological 
classes as a function of mass and distance from the cluster center. 


The completeness of this data sample is the same of our spectroscopic sample. This is defined by the 
ratio of the number of spectra yielding the redshift to the total number of galaxies in the photometric 
catalogue. |Cava et aTr| ( |2009| ) found that this sample is ~50% complete at V~ 18 mag. This corresponds to 
My ~ —19 mag and log(M*/MQ) ~ 10. The spectroscopic completeness rapidly decreases below V~ 19 
mag. The completeness turns out to be rather fiat with magnitude for most clusters, and it is essentially 
independent from the distance to the cluster center up to ~ 0.6i?200- 


Eig.[2 shows the distribution of the galaxy stellar masses, of the cluster-centric distance from the BCG 
normalized to R200, and of the absolute magnitudes for the different morphological types. Note that 
the histograms of the masses and luminosities remain fiat up to \og{M^ / Mq) ~ 9.5 and My ~ —18 
mag respectively. There are also no reasons to believe that the morphological fraction {i.e. the ratios 
among E-SO-Sp types) drastically change when the completeness decreases, and that the survey has 
systematically lost galaxies of a given morphological type in the explored interval of cluster-centric 


* The cluster membership was measured on the basis of the velocity dispersion of the galaxies in the clusters (see.jCava et al. 2009| 


^ i?200 is the radius delimiting a sphere with interior mean density 200 times the critical density of the Universe at that redshift (see for details 
[ 200 ^ 


Cava et al. 
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Figure 1. From top to bottom we plot here the distribution of the galaxy stellar masses, the cluster-centric distance in -R 200 units and the absolute magnitudes 
of the galaxies for the different morphological types: red lines mark the Elliptical galaxies, black lines the SO galaxies, and blue lines Spirals. The vertical 
dashed black lines mai'k the limit where the WINGS data sample is ~ 50% complete. 


distances 0 < Dbcg/R 200 < 0.6. We are therefore confident that up to these values of masses and 
luminosities the WINGS stacked sample is reasonably statistically robust. 

Fig. I^plots the distributions of R 200 , Cd (the central velocity dispersion of the cluster), log(Lx) (the 
X-ray luminosity of the cluster) and Md (the dynamical mass of the clusters) for the WINGS clusters 
that belong to the stacked sample. The values for the Coma clusters are indicated by the dashed lines as a 
reference. These data are useful to better figure out the observed range of properties of the stacked cluster 
sample. Clusters of quite different properties have been averaged together. 

2.3 Additional subsamples used in this work 


Three smaller subsamples have been used in this work for testing some of our hypotheses. The first 
one contains the measured Lick indexes for 2241 WINGS galaxies in common with our stacked sample. 
The Lick indexes were measured by Alexander Hansson in his PhD thesis following the prescriptions of 


Worthey et al. (19941. These data are still unpublished. The second subsample contains the bulge-to-disk 
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Figure 2. From top left to bottom right we plot here the observed distributions of -R 20 O’ ^ci for the WINGS clusters that belong to the 

stacked sample. The dashed line mark the values for the Coma cluster. 


decomposition of thousand of galaxies (but only 877 in common with our stacked sample) measured by 
Ruben Sanchez-Janssen in his PhD thesis. All galaxies were 2D decomposed with a Sersic bulge and 
an exponential disk component. These data are also still unpublished. The third subsample include 1579 


galaxies with redshift z < 0.1 extracted from the B-band PM2GC survey (see Calvi et al. [ 2011 2012 


for more details). These galaxies belong to low density environments and have been used to compare the 
properties of the different morphological classes in clusters and in the field. These galaxies have been 
analyzed adopting the same softwares (e.g. GASPHOT and MORPHOT) and procedures used in the data 
reduction of the WINGS data, so that there is a high level of homogeneity between the two databases. 
Clearly, all these subsamples are not complete from a statistical point of view, in particular for faint galaxies, 
but still they provide significant indications of the mean trends of the measured structural parameters. We 
will explicitly mention the use of these subsamples when they are used. 
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2.4 The high-z EDisCS stacked sample 


The properties of the EDisCS data survey ([White et al. [ [2005]) have been extensively analyzed in various 


papers (see e.g., Vulcani et al. 2011b). This multiwavelength photometric and spectroscopic survey 


includes several clusters at 0.4 < z < 1. The galaxy photometric redshifts were determined by using both 
optical and IR imaging ( jPello et ^T] |2009| ), the morphologies were derived by jPesai et al. [ ( [2007] ), and 
the photometric parameters (radii, total magnitudes, colors, etc.) were extracted using SExtractor by 


White et al. (2005). 


For the comparison with WINGS we used the mass-limited sample at M* = derived by 

Vulcani et al. | ( |20lTbl ), consisting of 489 galaxies, 156 of which are classified as Es, 64 as SOs, 252 as 
Spiral galaxies and 17 with unknown morphology. The values of R 200 ^nd the distance of the galaxies 
from the BCG were derived by jPoggianti et al. | ( |2006[[2()08| ). The values of the galaxy stellar masses (still 
unpublished) were kindly provided by B. Vulcani: they are calculated using the Kroupa ( Kroupa [ 2001) 
IMF in the mass range 0.1 — IOOMq following Bell & de Jong ( |2001 ). 


As for WINGS, the EDisCS data were also stacked together at increasing distance from the BCGs with 
the aim of providing a more robust statistical comparison between high-z and low-z clusters. 


3 OBSERVATIONAL ELEMENTS IN FAVOR OF THE TRANSFORMATION OF 
SPIRALS INTO SO S 


Here we present the main observational evidences in favor of the hypothesis that Spirals transform 
themselves into SO’s in clusters. 


Fig. shows the fraction of SO/E, SO/Sp and E/Sp for all the clusters observed up to z ~ 0.8, coming 
from the WINGS and EDisCS databases, as well as from literature data for the intermediate redshift 
(0.1 < 2 ; < 0.55) clusters studied by [Fasano et al. | ( |200^ and jPressler et al. | ( |1997| ). 


In their comparative analysis of the morphological fractions of clusters at low and high redshift Fasano et 


al. (2000) already claimed the existence of a general trend in the morphological fractions (%E’s, %S0’s, 


%Sp’s) and in the SO/E and SO/Sp ratios with redshift, clearly suggesting that the SO population tends to 
grow at the expense of Spirals, whereas the frequency of E’s remains almost constant. 


In Fig. [^the WINGS clusters are represented by one single point marking the average value found for 
our stacked clusters sample and its variance. There are several things that are worth noting in this figure: 1) 
the SO/E ratio (bottom panel) is almost constant with redshift with a significant variance; 2) the variance 
of the WINGS stacked sample is of the same order of the spread observed among clusters at all redshifts 
in all panels; 3) the SO/Sp ratio (mid panel) slightly increases at lower redshifts (z < 0.3). This fraction 
follows the SO/E ratio but with a smaller dispersion up to 2 ; ~ 0.15, where a strong peak is observed for 
two clusters A389 and A951; 4) the E/Sp fraction follows the same trend of the SO/Sp ratio, but with a 
bigger dispersion at all redshifts. Again a peak is observed at z ^ 0.15 for the same clusters; 5) the SO 
population dominate at z < 0.1. Depending on the cluster we can have a factor of ~ 3 — 4 more SO’s than 
Spirals; 6) at redshift around 0.1 and below the number of Spirals observed is very small. 


The spread observed in the morphological fractions can be explained remembering that: a) clusters have 
different structures and are more or less virialized/concentrated ( Oemler [ 1974| ); b) Spirals galaxies, that 
are the dominant population in the field, continuously enter into clusters from the cosmic-web filaments; c) 
there are possible errors in the classification of galaxies. 
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Figure 3. Trend of the morphological fractions (E/Sp, SO/Sp and SO/E) with redshift in various clusters. The WINGS stacked sample is shown with its Icr error 
bar. The plotted lines mark the evolution of the morphological fractions for a hypothetical cluster that at high redshift has the mean morphological fractions of 
the EDisCS sample. Spirals transform into SO’s and E’s at different rates per interval of redshift. Left panel: The solid lines mark the trend obtained for a rate of 
transformation of 6.25 Sp^SO’s per Gyr and 0 Sp^E’s per Gyr. The dashed lines instead are for a rate of 9 Sp^SO’s per Gyr and 0 Sp^E’s per Gyr. Right 
panel: The same plot, but with the following transformation rates: 5 Sp^SO’s per Gyr and 1.75 Sp^E’s per Gyr for the solid lines, and 7 Sp^SO’s per Gyr and 
1 Sp^E’s per Gyr for the dashed lines. 


loq(N 5 jj/N^^)=—0.83*loq(z} —0.74;rrns= 0.27:C.C =—0 82;Spear=-0,S9;Sign.= 1.53e-20 
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Figure 4. The log of the Nso/^Sp ratio versus the log of the redshift 2 . The solid line gives the best fit. Poissonian error bars mark the la uncertainty of the 
observed counts. 
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The general impression eoming from these data is that the effeets of the morphologieal transformation 
beeome visible at redshift lower than ~ 0.3, where a signifieant inerease of the E/Sp and SO/Sp ratios are 
observed. Despite the spread, the data seem to suggest an evolution with redshift of the morphologieal 
fraetions. If we trust in these data (i.e. that no systematie biases are present, e.g. in the morphologieal 
elassifieation and in the surveyed areaj^, we must eonelude from Fig.j^that there is an aetive meehanism 
transforming Spirals into SO objeets (and possibly into E’s). In faet, the ratios SO/Sp and E/Sp seem to 
inerease with approximately the same gradients, although with quite different dispersion. 

The lines shown in Fig. [^follow the hypothetieal evolution of a eluster eontaining 26 E’s, 15 SO’s and 68 
Spirals at redshift 0.9 (these values represent the average observed frequeneies of morphologieal types in 
EDisCS). They mark the results of a simple ealeulation in whieh Spirals are transformed instantaneously 
into SO’s and E’s with different rates per interval of redshift. In the left panel the solid lines mark the 
transformation rates of 6.25 Sp—>S0 per Gyr and 0 Sp— t-E per Gyr, while the dashed lines mark the rates of 
9 Sp—J-SO per Gyr and 0 Sp—;-E per Gyr. In the right panel we have instead for the solid lines 5.0 Sp—)-S0 
per Gyr and 1.75 Sp—;-E per Gyr, while for the dashed lines 7 Sp—)-S0 per Gyr and 2.5 Sp->E per Gyr. 
These rates were ehosen to fit the general trends observed and to demonstrate that we ean in prineiple follow 
the evolution of the morphologieal fraetion of any eluster aeross time. The observed trends demonstrate 
that an evolution of the morphologieal elasses is neeessary to explain the data. Our ealeulation however 
does not take into aeeount the infall, i.e. the faet that new Spirals enter the elusters aeross the eosmie time. 
Therefore, the rates ealeulated here should be eonsidered lower limits. 


Note that if we do not eonsider the transformation of Spirals into E’s, we eould not reproduee the mean 
trend of the SO/E fraetion: both the low (solid line) and high (dashed line) rates of Spiral transformation 
into SO’s ehosen in the left panel, in whieh the number of E’s remains eonstant, are not able to fit the mean 
WINGS point giving the SO/E fraetion. On the other hand, when the transformation into E’s is properly 
taken into aeeount (see the right panel), we ean reproduee both the mean trend of all morphologieal 
fraetions and the peaks in the E/Sp and SO/Sp ratio^as well as the nearly eonstant SO/E ratio observed 
for all elusters. In other word, or new E’s enter systematieally in the elusters or new E’s must be formed 
through major merging of Spirals. 


Fig. I^shows the log of the Nso/Nsp ratio versus the log of the redshift 2 ;. The least square fit is highly 
signifieant and demonstrate that the aetual data ean be interpreted in term of a temporal evolution of the 
morphologieal fraetion, in the sense that SO galaxies eontinuously inerease in number with respeet to 
Spirals. Transforming z in age of the Universe with a linear approximation in the redshift range 0 — 1, we 
see that the fraetion Nso/Nsp inereases by a faetor of ~ 2 every Gyr. 


In eonelusion, if our interpretation of Fig. and |^are eorreet. Spirals transform into SO’s and into E’s 
with a slightly different rate. Vuleani et al. | (201 lb), working on the same data for both surveys, arrived to 
the same eonelusion with their analysis of the morphologieal fraetions and mass funetions. 


Fig. [^provides us another hint in favor of the morphologieal evolution. It shows that SO’s and Spirals are 
not spatially distributed in the same way in elusters at low and high redshifts. The figure elearly demonstrates 
that in the WINGS elusters SO’s dominate the eentral region of the elusters (Dbcg/R 200 < 0.4), while 
Spirals in EDisCS are more frequent and abundant in the outer regions Dbcg/R 200 > 0.4. This faet. 


^ The cluster area observed is around 1 Mpc^ in almost all clusters. 

^ These peaks are not difficult to explain with a slightly higher frequency of transformation that can occur in some clusters in which the mechanisms of 
transformation are particularly efficient. 
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Figure 5. The left panels plot with a black line the difference between the number of SO’s and Spirals in the WINGS data sample for D^CG /-^200 <0.4 
(upper panel) and Dbcg/R 200 > 0.4 lower panel in each bin of galaxy mass. The right panels show the same difference for the EDisCS data sample. The red 
lines mark the same distributions for the difference between the number of SO’s and E’s. 


which reminds the morphology-density relation, is also eonsistent with the idea that Spirals transform into 
SO’s during their progressive infall toward the cluster central region. 

4 DISSECTING THE PROPERTIES OF THE STACKED CLUSTER SAMPLES 

In this section we analyze in detail the properties of the different morphological types as a function of 
mass and distance from the cluster center. The idea is that we can deduce from the observed differences of 
present day objects the mean path of the morphological transformation. In other words we can in principle 
establish which is the mean variation between two morphological classes. 

Fig. shows the results of our analysis for the WINGS stacked cluster sample. From top left to bottom 
right we have plotted here the following quantities derived for the different morphological types T in 
various bins of Dbcg/R 200 - the star formation rate (SFR), the log of the circularized effective radius 
(log(i?g[/cpc])), the mean effective surface brightness ({^)e[mag arcsec~^]), the log Sersic index (log(n)), 
the log of the luminosity weighted age (\og{wAge[yr])), the log total mass (log(M*/M0)), the B—V color, 
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Figure 6. From top left to bottom right we plot here versus the cluster-centric distance in R 200 units the following average quantities derived for the different 
morphological types in various bins of D^qq/R 2 Q 0 : the star formation rate (SFR), the log of the circularized effective radius ([og{R^[kpc]), the mean 
effective surface brightness ((/i}e), the log Sersic index (log(n)), the log of the luminosity weighted age ([og{wAge)), the log total mass (log(MQ), the B—V 
color, and the axis ratio (b/a). Red dots mark the Elliptical galaxies, black dots the SO galaxies, and blue dots Spirals. The horizontal dashed lines provide the 
average value of the morphological classes, while the colored strips mark the width of the error of the mean. The vertical dashed lines mark the limits where the 
stacked sample is ~ 50% complete. On the top of each diagram we report the mean difference between the SO and Spiral population. 


and the axis ratio (b/a). Red, blaek and blue eireles mark the mean value of eaeh quantity in eaeh bin 
of distanee respeetively for E’s, SO’s and Spirals. We deeided to plot with eolored strips the error of the 
mean a/'/N, where a is the standard deviation in the bin and N is the number of points, instead of 
the semi-interquartile in order to better highlight the mean differenees among the various morphologieal 
elasses. The width of the semi-interquartile is in faet ~ 3 times larger, so that in some eases the three 
morphologieal types share exaetly the same parameter. In eaeh diagram E’s, SO’s and Spirals are largely 
superposed, but the mean value for eaeh population is elearly distinguishable. 

Apart from the elear separation between SO’s and Spirals in eaeh plot, and the general similarity of SO’s 
with E’s rather than with Spirals (with the notable exeeption of the Sersie index n), there are two things 
that are worth noting in this figure. Eirst, SO’s masses are on average almost a faetor of ~ 2 larger than 
Spiral masses up to ~ 0.6i?20o; seeond, SO’s are ~ 25% smaller in their effeetive radius than Spirals within 
the same interval of eluster-eentrie distanee. 

Eig. |7] shows a different view of the same WINGS data. This time we averaged the same quantities per 
bin of galaxy mass. The same eolor eode marks again the three families of morphologieal types. Note that 
again SO’s are very similar to E’s in eaeh plot, and that at every galaxy mass E’s are systematieally eloser to 
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Figure 7. From top left to bottom right we plot here versus the log mass of the galaxies in solar units the following average quantities derived for the different 
morphological types in various bins of DbcgI^200- the star formation rate (SFR), the log of the circularized effective radius {log(R^[kpc]), the mean 
effective surface brightness in -R 200 units ((/i)e), the log Sersic index (log(n)), the log of the luminosity weighted age ([og{wAge)), the distance from the 
BCG in {Dbcg/^ 200 )^ the B—V color, and the axis ratio (b/a). Color codes as in Fig.^ The horizontal dashed lines provide the average value of the 
morphological classes, while the colored strips mark the width of the error of the mean. The vertical dashed lines mark the limits where the stacked sample is 
~ 50% complete. On the top of each diagram we report the mean difference between the SO and Spiral population. 


the cluster center than SO’s and Spirals. This is a manifestation of the morphology-density relation (Dressier 
[ 1980 ). Furthermore, SO’s are at each mass on average ~ 25% smaller than Spirals, as seen before. Note 
that if we restrict the sample to a given threshold, e.g. in total mass, we observe that the average values 
shown in each plot slightly vary, but the relative differences among the various morphological classes 
poorly depend on the adopted cut-off level. This suggests that the result is statistically quite robust and not 
strongly related to the completeness of the data sample. 


To test the statistical robustness of this result we carried out a bootstrap resampling ( Efron & Tibshirani 
d 1993[ ) aiming at measuring the significance of the difference in the mean mass values for galaxies of 
morphological type SO and Sp. The bootstrap analysis confirms that the difference between average M* 
of SO and Sp is highly significant. Only in 1% of 10“^ bootstrapped samples the difference between SO 
and Sp M* averages is significant at a confidence level lower than 2 ct. The two samples are different also 
according to a Kolmogorov-Smirov test carried out on each bootstrapped sample: again only 1% of 
the bootstrapped samples show Kolmogorov-Smirov’s D that implies a significance lower than 2cr. This 
result is valid even if a restriction to the cluster-centric distance r < 0.2i?200 is applied. It is stronger if 
r < 0.6i?200 because of a larger original sample size. It is also valid if the initial mass of the SSP at age 


zero, or the mass in a nuclear- burning phase of (as defined following Fritz et al. 20111 are considered. 
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An additional confirmation is offered by the mass computed from surface photometry Mp^, following |Bell| 
& de Jong ( |2001 ). In this case there are significant differences with M*, which for the full WINGS SO 
sample follow a law 5 log M = log M* — log Mph ~ —4.03 + 0.688 log M* — (0.0296 log M*)^, yielding a 
6 log M ^ —0.25 at log M* 9. A similar effect is also seen for E and Sp galaxies, and remains present if 
restrictions to r < O.2i?200 and r < 0.6i?200 are applied. Nonetheless, a bootstrap applied to Mph confirms 
the same difference and the levels of significance as for M*. 


The differential distribution of the masses for the three morphological types in the different bins of 
cluster-centric distance clearly indicate that the higher masses (> lO^^M©) are almost exclusively SO’s 
and E’s in the 0 — 0.2 Dbcg/R 200 bin. In the 0 — 0.2 interval of Dbcg/R 200 we count 259 E’s, 395 SO’s 
and 135 spirals, so the statistics here is quite robust. At larger distances the number of galaxies of different 
morphologies become approximately equal, but the number of E’s rapidly decreases. Then, in the last bin 
at 0.8 — 1.0 Dbgg/R 200 , the lack of spatial completeness of the spectroscopic survey starts to be severe, 
so we cannot trust very much on the observed distribution. 

Eig.j^and Eig.|^show how the fractions Nso/Nsp and Nso/Ne vary as a function of different variables 
in various bins of cluster-centric distance. We see that there are ~ 50 times more SO’s than Spirals with 
mass > IO^^Mq and up to ~ 20 times more SO’s than E’s of low masses in the first bin. The ratio 
Nsq/Nb increases at low masses, but this result is more uncertain being our sample less complete below 
log M* /Mq < 9.5. In Eig. [^the distribution of the effective radii Re looks bimodal for the SOs/Sp ratio, 
indicating that many of these galaxies are massive, but likely of small dimension, i.e. quite dense objects, 
as revealed by their higher surface brightness. They are probably bulge-dominated objects. It is interesting 
to note the distributions of the Sersic index, of the luminosity weighted age, of the B—V color, and of 
the mean effective surface brightness. All these structural and photometric parameters indicate that in 
the 0 — 0.2 interval of Dbgg/ R 200 there is a well defined population of SO’s with n values larger than 
those of Spirals and lower than those of E’s, preferentially older than Spirals, but younger than E’s, redder 
than Spirals but bluer than E’s, with a higher surface brightness with respect to Spirals, but similar to E’s 
in all the bins. We also observe that the SO’s with small radii have also low n values, young ages, blue 
colors and SERs similar to Spirals. Going at larger cluster-centric distances we instead see a clear trend 
toward Nso/Nsp ~ 1 and Nso/Nb ^ I in almost all diagrams with the obvious exception of the surface 
brightness and the axis ratios. The general impression is that all the galaxies closer to the cluster center 
are systematically different from the others. This fit with the idea of a progressive build up of the cluster, 
where the most processed objects (red and dead) are the first entered in the cluster core, while those in the 
outer parts are still in a phase of rapid transformation showing traces of residual activity. 


Looking at the above figures we must conclude that the effeets of galaxy transformation should determine 
the following changes in the past Spirals (now SO’s) with respect to present day Spirals: 1) the mean age, 
mass and color, 2) the concentration of the light profile so thaf smaller effective radii and larger surface 
brightness are measured, 3) the shape of the light profile so that the Sersic index increases, 4) progressive 
cut-off of the SER. We will discuss these possibilities in Sec.|^ 

Eig.[T^and Eig. [^present the behavior of the EDisCS galaxies in various bins of cluster-centric distance 
and galaxy mass as we did in Eig. |^and Eig. [^respectively for the WINGS data. These plots have been 
obtained in the same way by stacking the EDisCS data as we did for the WINGS sample. Here again the 
underlying hypothesis is that the properties of the clusters at high redshifts are quite homogeneous and 
the effects of galaxy evolution not too heavy to change systematically the properties of galaxies in the 
redshift interval considered. Unfortunately the EDisCS data have not the same degree of accuracy and 
quality, therefore they are used here only for a qualitative comparison with WINGS. 
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Figure 8. From left to right: the log of the total galaxy mass (log(MQ), the circularized effective radius {\og{R^[kpc]), the log Sersic index (log(n)), and the 
log of the luminosity weighted age (log(u;A^e)). The y-axis shows the ratios between the number of SO’s and Spirals (black line) and the number of SO’s and 
E’s (red line) in five different ranges of distances Dqcg /^200 (0-0.2, 0.2-0-4, 0.4-0.6, 0.6-0.8, and 0.8-1.0). 


With respect to WINGS (limited to the same interval of galaxy masses) we note the following behavior at 
each cluster-centric distance : 1) the masses of SO’s and Spirals are nearly equivalent (the mean difference 
being lower than 0.1 dex for both samples); 2) Spirals are still systematically larger than SO’s in both 
samples (5*0 — Sp = —0.08 in EDisCS and SO — Sp = —0.14 in WINGS). This result however is biased 
by the fact that the EDisCS effective radii are not deconvolved for the seeing point spread function; 3) in 
WINGS there is a mean systematic difference in the mean surface brightness of SO’s and Spirals of ~ 0.58, 
while it is only ~ 0.1 in EDisCS. Again this result is a direct consequence of the measured effective radii; 
4) Spirals are systematically brighter than SO’s in EDisCS of ~ 0.5 mag, while the difference decreases to 
~ 0.2 in WINGS; 5) the difference in the mean B — V color is nearly equivalent (~ 0.1 mag), while the 
EDisCS sample is systematically bluer as expected on the basis of their redshift, 6) there is a substantial 
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Figure 9. From left to right: the B—V color, the star formation rate (SFR), the mean effective surface brightness ((/i)e), and the b/a ratio. As in the previous 
figure the y-axis gives the ratio between the number of SO’s and Spirals (black line), and SO’s and E’s (red line). 


difference in the mean flattening of Spirals and SO’s that we believe could be connected to the difficulty of 
detecting and classifying nearly edge-on objects at high redshift. 

Notably in each bin of galaxy mass the cluster-centric distance of Spirals with respect to SO’s and 
E’s is much higher (50 — Sp = —0.23 in EDisCS and SO — Sp = —0.07 in WINGS). Spirals are also 
systematically brighter than SO’s up to ~ lO^^M© in EDisCS (50 — Sp = 0.51 vs 0.16 in WINGS). 

The general impression from these figures is that even at high redshifts the population of SO galaxies is 
already quite different from that of Spirals. Once the effects of passive evolution are taken into account, we 
can say that the Spirals and SO’s at high redshift are quite similar to those observed today, apart from their 
relative frequency. This means that, either it exists a population of primordial SO’s, or the transformation of 
Spirals into SO’s was active since the first epoch of cluster formation. Furthermore, if a transformation has 
occurred it should have been quite rapid, probably not longer than ~ 1 — 3 Gyr, because we are not able to 
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Figure 10. The trend of the measured parameters for the EDisCS sample as a function of the cluster-centric distance. 

detect a class of objects with mixed properties between Spirals and SO’s in both surveys. See discussion in 
Sec. HI 


5 OBSERVATIONAL EVIDENCES SUPPORTING THE GAS STRIPPING BY A HOT 
MEDIUM 


Ram pressure stripping is one of the physical mechanisms invoked to explain the fast transformation of 
Spirals into red passive galaxies similar to many SO’s. Direct evidences that this phenomenon is currently 
active in the WINGS clusters have been presented by Poggianti et al. ( 2015| ). They have identified 
241 possible Jellyfish galaxies that exhibit material outside the main body of the galaxies, suggesting a 
phenomenon of gas stripping. 
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Figure 11. The same trends as a function of galaxy masses. 


Another observational faet suggesting that Spirals are affected by hot IGM interaction comes from Fig.[T^ 
that shows the relationship between the X-ray luminosities (in log units) of the WINGS clusters and their 
morphological fractions SO/Sp and E/Sp. Here we have taken into account only those WINGS clusters 
with a robust number of detected morphological types (38 clusters). We have excluded those clusters with 
less robust statistics in which the number of Spirals members is < 2. The figure clearly indicates that a 
significant degree of correlation is observed in particular for the SO/Sp fraction. The most luminous X-ray 
clusters have preferentially a larger SO/Sp (and possibly E/Sp) fraction. This is a clear hint that in the hotter 
IGM the stripping of gas and consequently the morphological transformation of Spirals is more efficient. 
The gas is stripped by the hot IGM when Spirals enter the cluster environment and start orbiting around the 
cluster center. 
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Figure 12. The panel shows the X-ray luminosity of the WINGS clusters plotted against the morphological fractions SO/Sp (upper panel) and E/Sp (lower 
panel). Poissonian error bars have been assumed. 


The dependence of the morphological fractions on the X-ray luminosity of the clusters, and consequently 
on the temperature of the intra-cluster medium, led us to suspect that the mechanism of evaporation 
(introduced by|Cowie & Songai'Ia~l|1977[) might play a role as significant as that of ram pressure. 


In general, if Spirals transform into SO’s via gas stripping and quenching of SF, one might suspect that 
the mass distribution of present day SO’s should not be too different from the mass distribution of early 
Spirals. Loosing their gas Spirals do not change very much their total mass (< 10% variations are expected). 
This is indeed observed in Fig.[^ This figure plots the mass distribution for the two samples of WINGS 
and EDisCS for galaxy stellar masses above 10 ^^^'^Mq. Note the strong similarity between the WINGS 


SO’s and EDisCS Spiral population. This similarity is also very clear in Eig. 7 of Vulcani et al. (20111 
showing the mass function of the two morphological classes in WINGS and EDisCS. Certainly this is not a 
conclusive argument, but it contributes to the general picture in favor of the galaxy transformation. 


A further evidence that galaxies are stripped of their gas, particularly in the central region of the clusters, 
comes from Eig. [^that shows the behavior of the Lick index for the WINGS galaxies measured by 
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Figure 13. From top to bottom the mass distribution of the WINGS SO galaxies compared with that of Spirals, E’s and SO’s from the EDisCS galaxies. 


Alexander Hannson in his PhD thesis. The figure elearly indicate that the galaxies in the outer region of the 
clusters still present this line in emission, while the galaxies in the vicinity of the center are systematically 
deprived of their gas. A more detailed analysis of the spectral emission properties of the WINGS galaxies 
will be presented by Marziani et al. (in preparation). 

Are these hints sufficient to claim that the gas stripping is the only mechanism producing the 
transformation of Spirals into SO’s? Certainly not. In the next section we will see that do exist several 
suspects that merging events and gravitational interactions are active somewhere in the clusters and can 
contribute to the formation of SO galaxies. 


6 OBSERVATIONAL EVIDENCES SUPPORTING THE MERGING HYPOTHESIS 


As mentioned in the Introduction, the possibility that SO galaxies are produced via the merging of two 
spirals with unequal mass has been discussed since the late 1990s. Bekki | ( [1998 l found that simulations 
of the merging of two bulgeless galaxies (with the gas component modeled) and a mass ratio 0.3 could 
reproduce the photometric and structural properties of a typical SO galaxy (red color, absence of a stellar 
thin disk). [Aguerri et al. (2001) discussed a similar scenario for SO formation, and found a steepening 
of the photometric profiles toward the early-types due to satellite accretion. After a minor merger, bulges 
grow and the Sersic index increases. Eliche-Moral et al. (2013) found that intermediate mergers (i.e., 
events with mass ratios between 1:4 and 1:7) can produce remnants that photometrically and kinematically 
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Figure 14. The Lick index for galaxies closer to the center of the clusters and for more distant objects. 


resemble SO galaxies, eonsistent with the bulge-to-disk coupling observed in real SO’s. A significant disk 
heating is expected from minor mergers (|Qu et al.' 20111. 


These models provide support to the idea that minor/medium mergers lead to an increase of the Sersic 
index value, if spheroids or B+D systems are involved. The simulations of Hilz et al. ( 2013[ ) show 
that unequal mergers (m 2 < 1/5) produces a steep increase in the Sersic index in massive spheroid 
(log M > 11). Repeated minor merging can even lead to the build-up of elliptical galaxies without resorting 
to major mergers ( [Boumaud et al. 20051. [Querejeta et al. (20151 computed photometric profiles from the 
dissipative N-body binary merger simulations from the GalMer database, and found that minor merger 
remnants can have low Sersic indices 1 < n < 2, consistent with the presence of pseudo-bulges in 
observed SOs. In general therefore N-body simulations predict that the Sersic index of the galaxies evolve, 
proportionally increasing with the dimension of the new formed bulge component. 

From an observational point of view we see evidences of ongoing tidal interactions in the WINGS 


galaxies (see e.g. Poggianti et al. , 2015). This means that gravitational interactions and merging can still 


be active today in particular in the outskirts of clusters. Unfortunately the analysis of present day SO’s does 
not easily reveal whether they are the end product of past merging events. Direct proof of past accretion 
events could come from the analysis of the frequency of GCs or PN, or from a detailed analysis of the 
stellar populations and SFH. This is quite difficult already for very nearby objects, but almost impossible 
for far SO’s that can be studied only in their integrated properties. 
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Figure 15. The bulge-to-total luminosity ratio B/T, the effective radius of the bulge, and the Sersic index of the WINGS SO galaxies are plotted here versus 
the K band absolute magnitude. The red and green dots mark the bulge-to-disc decomposition of the merging products coming from N-body simulations of 
IQuerejeta et al. |{2Q15| and |Tapia et al. |(2Q14) respectively. 


The only possible approaeh with our data is that of statistieally distinguish the mean properties of the 
different morphologieal elasses, and to eheek whether the observed differenees might be explained by 
theoretieal models. Fig. 15 eompares the B/T ratio and the Sersie index derived from N-body simulations 
of major and minor merging events among binary galaxies with the eorresponding values measured for the 
WINGS data set by Sanehez-Jannsen in his PhD thesis. The figure suggests that the formation of eluster 
SO galaxies by merging eould be one element to eonsider, at least regarding its brightest representatives. 
However, it is not elear yet whether smaller eluster SOs ean be the end produet of merging events. 


The problem is that the high veloeity dispersion in elusters does not favor the oeeurrenee of merging. We 
therefore suspeet that most of the merging events happen when galaxies reside in the outskirts of elusters or 
when they are still members of small groups that are infalling into the elusters. This means that a substantial 
pre-proeessing eould be experieneed by Spirals before being virialized into the eluster environment. 

A seeond hint in favor of the merging hypothesis eomes from Fig.[^ where we have seen that the fraetion 
SO/E at low redshift eannot be explained without a small inerease of the E types, and this might be explained 
with major/minor merging events. 


Therefore, prompted by the idea that merging events are neeessary for explaining the aetual strueture of 
SO galaxies and the frequeney of the morphologieal fraetions in nearby elusters, we deeided to eheek if 
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Figure 16. The Sersic index as a function of galaxy stellar mass for the three morphological types in two different regimes of cluster-centric distances. Black 
dots mark the values for Dbcg/^200 ^ 0.3 while red dots those for Dbcg/R 200 > 0.3. 




Figure 17. Left panel: The B/T ratios of the various morphological types in the two bin of cluster-centric distance. Right panel: The Sersic index of the bulge 
components of the galaxies. 


the WINGS data in general, and the Sersie index of SO’s and Spirals in partieular, are eonsistent with the 
merging hypothesis. If sueh hypothesis is eorrect, we expeet to see an evolution of the Sersie index, i.e. 
that n inereases as a eonsequenee of multiple merging events experieneed by Spirals. 
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In Fig. 1^ and 1^ we have already shown that the Sersie inde?0of the various morphologieal elasses are 
quite distinet both for various masses and for distanees from the eluster eenter. The two figures showed 
that for a relatively broad mass range (9.5 < log{M^ /Mq) < 11) the mean Sersie index of SO’s and today 
Spirals has a systematie differenee of A log(n) ~ 0.2. In other words, in the same mass range SO’s and 
Spirals have different light profiles. We ean distinguish three broad mass domains: small stellar mass 
log(M*/MQ) < 9.5: the systematie differenee between SO’s and Spirals is small. In this mass domain, we 
have the least massive systems with the smallest bulges, and n is eonsistent with the presenee of pseudo 
bulges, i.e. of struetures that might originate from the seeular evolution of disks. In the intermediate mass 
range, we observe instead systematieally higher n for SO’s and E’s with respeet to Spirals. For masses 
\og{M^/ Mq) > 11 the n value for Spirals beeomes not reliable beeause there are too few objeets of this 
mass. The general trend is that as far as the spheroidal eomponent of the galaxies inereases, the value of n 
inereases. 


Trusting in N-body simulations we expeet that the larger Sersie index of SO’s and E’s at eaeh mass 
eould be originated through a series of major/minor merging events experieneed by a Spiral galaxy. Ram 
pressure and evaporation eannot in faet drastieally ehange the Sersie index of a galaxy. A redistribution or 
an aeeretion of new stars is required to produee the observed variations, in partieular of the bulge-to-disk 
ratios. The transforming galaxies should therefore experienee sueh merging events or should produee new 
stars in their bulges. 


Vuleani et al. | (2011b) eame to a similar eonelusion in their analysis of the mass funetions and 


morphologieal mix of the WINGS and EDisCS data. They stressed that it is impossible to reeoneile 
the mass funetion observed at low redshift only with the evolution of the morphologieal fraetions. Mass 
aeeretion events and new infalling galaxies are required to reproduee the observed evolution. They believe 
that new SE events should oeeur in the transforming galaxies. 


A robust hint that merging events and gravitational interaetions are responsible of the observed n — M* 
relation would be the deteetion of a systematie ehange of n in galaxies that potentially might be subjeet 
to several eneounters. We have seen before that the eenter of elusters, where the density of galaxies and 
the veloeity dispersion is high, host the more massive objeets (a faetor of ~ 2 larger in mass) as well as 
the smaller in size (~ 25% smaller, with the notable exeeption of the BCGs). This differenee eannot be 
explained in the eontext of ram pressure stripping and evaporation by a hot IGM, but fits with the idea that 
a series of merging events and strong gravitational interaetions have produeed the big SO’s that are now 
seen in the eluster eenter. 


In order to test sueh hypothesis, we examined the behavior of the Sersie index of the three morphologieal 
elasses in Eig. 16 Here for eaeh elass we plotted the mean galaxy mass and Serise index observed in two 
intervals of eluster-eentrie distanee Dbcg/R 200 < 0.3 (blaek dots and gray strips) and Dbcg/R 200 > 0.3 
(red dots and blue strips). The figure shows that there is only a marginal evidenee of an inerease of n for 
the SO’s in the eluster eentral region. The probability that the observed means are different by ehanee is 
~ 10% and the observed variation is very small. On the other hand the figure shows that inner and outer 
E’s have signifieantly different values of n. In partieular the inner E’s have systematieally lower n than 
outer objeets. The statistieal signifieanee that the two samples have different averages by ehanee is now 
~ 1%. This means that E’s in the eentral region of the elusters are likely subjeet to tidal stripping (galaxy 


^ The Sersie index was measured by GASPHOT (see, |D’Onofrio et al.; |2Q14| on the major and minor growth curves of the galaxies without distinguishing the 
bulge and disk components. 
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harassment), i.e. the most external stars are lost producing less extended envelopes and consequently lower 
measured n. Spirals on the other hand remain the same at every distance inside the clusters. 


One might speculate that the variations of n could not be associated with a real change of the bulge 
structure since GASPHOT is a software that provides the Sersic index for the whole galaxy and not for the 
bulge component. For this reason we used again the subsample of galaxies of Sanchez-Jannsen, which 
contains the bulge-to-disk decomposition of our galaxies. Fig. [TT] confirms that: 1) there is a systematic 
variation between the Sersic index of the bulges of the various morphological types; 2) the mean difference 
between the B/T ratios of SO’s and Spirals increases in the inner region of the clusters; 3) the mean 
difference in the Sersic index also increases. 


In conclusion we can say that we have only a marginal indication that the SO’s bulges are systematically 
bigger in the central region of the clusters as a consequence of major/minor merging events, but in general 
the data do not support the idea that merging events are more frequent in the inner region of clusters with 
respect to the outskirts. The question is therefore when and where the transforming Spirals acquire the new 
stars that complete the transition to the SO class systematically increasing their bulge structures. 


Up to now several works have already claimed that the main quenching and morphological transformation 
do not occur primarily in the cores of clusters, at least not at z < 1. The decrease in SF is already visible at 
several virial radii ([Lewis et al. j |2002t [Gomez et al. i|2003|), and the SO fraction increases since z ^ 0.5 


(Dressier et al. 1997; Fasano et al. 2000) particularly in less-massive clusters (Poggianti et al. 2009 


Just et al. 2010). The evolution driven by the environment seems to occur at intermediate densities and 


outside the cluster virial radius and such effects are observed in simulations out to as many as five virial 
radii [Bahe et al. [ (2013). This open the possibility that pre-processing mechanisms affect the infalling 
galaxy population. 


The emerging scenario from the above discussion is that Spirals start to quench their SF when they 
enter in the hot IGM, but at the same time when they are still in the cluster outskirts they progressively 
accrete/merge new material changing the structure of their bulge and disk (thick disks?). One possibility 
is that the Spirals which end up in today SO’s, were members of small groups infalling the clusters and 
subsequently merged together. 


In other words from the present sample we can only derive indirect hints that merging events have 
contributed to the formation of today SO’s. This il likely due to the fact that the present data cover only 
the inner virialized part of the WINGS clusters. Probably the ongoing observations of the outskirts of the 
WINGS clusters will provide soon much solid conclusions about this. 


An intriguing fact that is difficult to explain with the ram pressure is that the clusters with a small value 
of the Nso/Nsp ratio (< 1) seem to contain in their central region systematically more massive SO’s than 
clusters with a higher value of this ratio (see left panel of Fig. [T8| ). We have verified that there is not a big 
difference in the mean luminosity weighted age between these two types of clusters (~ 0.1 dex), so we 
cannot attribute the effect to a systematic variation in age. 


However, Fig. 19 suggests a partially significant correlation between the log of the Nso/Nsp ratio 
and the central velocity dispersion of the clusters. A bisector least square fit gives: \og{Nso/Nsp) = 
1.64 log((Tc^) — 4.42 with a rms = 0.2 c.c. = 0.42 and a significance of 1.79 10“^). A weighted least 
square fit provides a shallower slope: \og{Nso/Nsp) = 0.87\og{aci) — 2.26. On average the fraction 
Nso/Nsp increase by a factor 2 — 3 going from ~ 500 to ~ 1100. This implies that the more 
massive galaxies reside in the clusters with the lowest velocity dispersion, in agreement with the idea that 
merging events are more frequent because the velocity dispersion in the clusters is low. 
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Figure 18. Left panel: The masses of the three morphological types in clusters with a high and small ratio of SO’s to Spirals. Right panel: The distribution of 
the masses of the three morphological types in clusters with a different fraction of SO/Sp. 


noted that the Nso/Nsp ratio is lower in more eoneentrated elusters where the 
number of E’s is higher. The suspeet is then that interaetions and merging have played a key role in 
determining the high mass of these galaxies. Probably in these elusters the merging aetivity was higher 
and the Nso/Nsp ratio kept small either by the lower temperature of the X-ray gas (remember that the 
Nso/Nsp fraetion depends also on Lx) and by a eontinuous aeeretion of new Spirals from the eosmie web. 

The right panel of the same hgure shows that there are no systematie differenees in the mass distributions 
of the various morphologieal types in elusters with Nso/Nsp < 1 and Nso/Nsp >2.5 suggesting that the 
observed behavior is not the result of a seleetion bias, but is probably an intrinsie eharaeteristie of galaxies 
of these two types of elusters. 

At the end of this Seetion we should say that we have aeeumulated several hints, although none strietly 
eompelling, that merging events and gravitational interaetions (harassment?) have had a role in galaxy 
elusters (in partieular in those with low veloeity dispersions). 


Fasano et al. (2000 


7 FIELD AND CLUSTER SO’S 


One of the most relevant things of the previous analysis is that E’s in the eentral region are systematieally 
deprived of their outer stars with respeet to more distant objeets. This tell us that gravitational stripping and 
galaxy harassment are playing an aetive role in the high density regions of elusters. Galaxy harassment may 
be assoeiated with the origin of a diffuse stellar eomponent that yields the so-ealled intra-eluster optieal 
light ([Zwieky |(|1952[)), if this diffuse eomponent is due to stars tidally stripped from galaxies as they move 


within the eluster (e.g., Gregg & West 1998). 


To test this effeet we eompare now the properties of the three morphologieal elasses for objeets belonging 


to the general held and to elusters. We have therefore used the PM2GC data sample deseribed in See. 2.3 


Fig. shows the eomparison of the Sersie index and the effeetive radius of the three morphologieal types 
for galaxies in elusters (blaek dots and gray strips) and in the held (red dots and blue strips). The left panel 
tells us that Spirals in elusters and in the held have exaetly the same n. On the other hand, SO’s and E’s 
in elusters appear systematieally with lower n values at eaeh mass. We have eheeked with a Weleh’s test 
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Figure 19. The panel show the log of the Ngo /Ngp fraction against the central velocity dispersion of the clusters. The solid line is the bisector least square fit, 
while the dashed line is the weighted least square fit. 


that this difference is significative from a statistical point of view and cannot originate from the different 
wavebands (B and V) of the two datasets. The mean difference that can be attributed to the waveband effect 
is of the order of ~ —0.05 for A log(i?e[^ — B]) and ~ —0.2 for An[V — B]. The right panel confirms 
that the higher Sersic index of galaxies in the field is associated to a systematic variation of the circularized 
effective radius. Field galaxies are systematically larger than cluster objects, suggesting that the main 
visible action of the cluster environment is the stripping of the outer regions of galaxies. 

The main indication is that the mean E’s and SO’s of intermediate mass in the field have larger n 
(~ 10%) and larger Rg (~ 20%). One possibility is that during the infall into the cluster potential, or 
as a consequence of repeated high speed encounters, galaxies loose their outer regions. This produces 
systematically smaller radii and also smaller Sersic indexes for the galaxies in clusters. Moreover, the 
similarity of the n — M* relation of field and cluster galaxies suggests that the mechanisms that have 
produced the various morphological classes should be very similar in both environments. 

If ram pressure and evaporation were not active in the field, the only possible hypotheses are that field 
SO’s originated from Spirals that were members of small groups that are now dissolved, or in binary 
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Figure 20. Left panel: The Sersic index of E’s (upper panel), SO’s (middle panel) and Spirals (lower panel) of galaxies in clusters (back dots and gray strips) 
and in the field (red dots and blue strips). Right panel: The effective radius in log units. Symbols as in the left panel 


associations. Part of these objects could also be byproduct of a secular evolution of the disk. In other words 
the merging hypothesis is still one the most attractive even for field SO’s. 


8 CONCLUSIONS AND DISCUSSION 

In this work we have obtained the following results: 

1. the fractions SO/Sp and E/Sp derived for clusters at different redshifts smoothly increase from z ~ 0.8 
to the present, while the fraction SO/E remains almost constant (or is only slightly increasing) although 
with a significative variance. In particular the fraction SO/Sp is nearly double every Gyr; 

2. starting from the mean number of morphological types observed in EDisCS and ending to those 
measured in WINGS, we calculated a mean rate of transformation of Spirals into SO’s of ~ 5 objects 
per Gyr (not considering the infall of new Spiral objects in the clusters), while the mean transformation 
of Spirals into E’s is at least ~ 2 objects per Gyr; 

3. the SO population in clusters appears to grow at the expenses of the Spiral one by a factor of ~ 3 — 4 
since 2 ; ~ 1. At the same time the number of E’s could not remain constant; 

4. SO’s are more frequent in the central region of nearby clusters, while Spirals are preferentially found in 
the outer regions of high redshift clusters; 

5. the SO family at z ~ 0 is on average quite different from that of present day Spirals: in general at a 
given distance from the cluster center they are more massive (up to a factor of 2), 25% smaller in size, 
with an higher Sersic index and effective surface brightness, and are redder in B—V colors. On the 
other hand, at a given mass they are closer to the center, systematically smaller in size, with higher 
surface brightness, age, color and Sersic index, but lower SER; 

6. the more massive SO’s in WINGS are clustered around the center of the clusters; 

7. the SO’s in EDisCS have the same structural characteristics of present day SO’s onee a passive luminosity 
evolution is taken into account. As for the WINGS clusters E’s and SO’s are more frequent in the 
central part of the clusters, but the proportion of the morphological types is quite different; 
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8. the mass distribution of present day SO’s and high redshift Spirals is very similar (although in the 
limited range of galaxy masses M* > 10 

9. the morphologieal fraetions SO/Sp and E/Sp in the loeal WINGS sample are mildly eorrelated with 
the X-ray luminosity of the elusters and with the eluster veloeity dispersion. The objeets eloser to the 
eluster eenter are systematieally deprived of their gas. The more massive SO’s reside in the elusters 
with the smaller veloeity dispersion; 

10. the Sersie index of Spirals, SO’s and E’s inerease with the mass of the galaxies and in partieular with 
the bulge eomponent. At eaeh mass the value of n is large for E’s, intermediate for SO’s and small for 
Spirals; 

11. Pairwise merger simulations give some support to the idea that the most massive eluster SOs may arise 
by merging; 

12. there are no elear indieations of an inerease of the Sersie index for galaxies eloser to the eluster eenter 
with respeet to more distant objeets, with the possible exeeption of eentral E’s (BCGs exeluded) that 
appear to have systematieally lower n of far E’s; 

13. the eomparison between eluster and field galaxies tells us that galaxies in elusters are systematieally 
smaller in effeetive radius and have lower n values. 

The first four points of this list strongly suggest that Spirals transform into SO’s and E’s in elusters. 

The fifth and seventh points eonfirm that SO’s and Spirals are very different in many observable parameters 
at all redshifts. If we aeeept the idea that Spirals transform into SO’s, we should reeognize that the 
morphologieal transformation is radieal. The speed of the transformation should be quite fast, sinee the 
various surveys at different redshifts have elearly identified the three morphologieal elasses with almost no 
objeets with eonfused or mixed morphology. 

The eighth point suggests that the mass distribution of present day SO’s and Spiral is eonsistent (at least 
for the large masses) with the hypothesis of transformation through gas stripping and quenehing of SE. 
However, within this seenario it is quite diffieult to aeeount for the large observed differenee in the Sersie 
index among the progenitor Spiral and the final SO galaxy. 

The ninth point tell us that the hot IGM eertainly play a key role in the transformation beeause the elusters 
having a bright X-ray luminosity eontain proportionally more E’s and SO’s with respeet to Spirals. However, 
the same point tell us that the elusters with the smaller veloeity dispersion eontains the biggest SO’s in 
agreement with the idea that merger events have produeed these galaxies. 

The aetivity of merging and gravitational interaetions is suggested by points 3, 10 and 13. The main 
hints that SO’s eould originate from merging events (if we aeeept the idea that Spirals must transform), rest 
on: a) the faet that up to Dbcg/R 200 < 0.7 at eaeh distanee form the eluster eenter the mean Spirals are 
systematieally less massive than the mean SO’s and E’s; b) in the same interval the mean E’s are more 
massive than the mean SO’s and the most massive are progressively segregated in the eentral region; e) the 
Sersie index is predieted to evolve after merging and is observed to inerease systematieally with the mass 
of the spheroidal eomponent (and marginally deerease with the distanee from the eluster eenter). 

The twelfth point indieates that the merging aetivity is not eonfined to the densest regions of the elusters, 
but ean oeeur in every plaee and time inside them, in partieular in the outer eluster regions, where the 
veloeity dispersion slow down. We ean therefore speeulate that when small groups fell into the eluster 
potential they merged together forming the new SO’s and E’s we see today. 
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The thirteenth point suggests that the main action of the cluster environment is that of stripping stars from 
the outer parts of the galaxies. All the morphological types in clusters appear in fact to be systematically 
smaller in size and with lower Sersic index with respect to field objects. This effect of harassment is not in 
contradiction with the idea that merging is the main actor of the Spiral transformation. 


The merging hypothesis is currently supported by recent observations of massive SO’s located in groups 
and clusters, usually exhibiting traces of past merger events ([Rudick et al. [|2009t|Janowiecki et al. j|2010 


Mihos et al. [ 2013|). Furthermore, the high resolution analysis of the AtlasSD group (Cappellari et al. 


2011 2013 Emsellem et al. 20111 has found that 87% of the galaxies classified as SO’s or E’s at very low 


redshifts are fast rotating systems, indicating that they have rotationally supported disks. Slow rotators are 
the most massive galaxies, often with boxy isophotes (70% of them), and depletion of mass in the central 
galaxy regions (core). East rotators are less massive, typically have disky isophotes, and centrally peaked 
surface brightness profiles (cuspy). East rotators are presumably stripped Spirals with some mass added by 
minor mergers, whereas the slow rotators are formed in more violent processes so that the bulges formed 
first. 


This scenario is also consistent with other well known observations: 1) stellar and gaseous disks, and 
dust lanes have been frequently found in E’s and SO’ bulges; 2) shells and ripples are often seen around 
these galaxies when deep exposures are taken; 3) lenticular galaxies surrounded by polar rings exist; 4) 
early-type galaxies could evolve in size by multiple minor merging events (see g.g., |Lopez-Sanjuan et al. [ 


20121 . 


All these facts together point toward a strong role of merging in re-structuring the morphology of late-type 


galaxies. This is supported by modem high resolution simulations of merging events (e.g. Cox et al. 2006 


Robertson et al. [ 20061 between two Spiral galaxies that have produced E’s and SO’s with some rotational 
support with the right coupled bulge-disk structures in time scale less than ~ 3 Gyr ( Querejeta et al. [ 20151. 
A major merging could completely transform the structure of a Spiral galaxy, varying its bulge-to-disk 
ratio, its Sersic index, the mean surface brightness and the effective radius (see g.g., [Querejeta et ah] 20151. 
Minor mergers can also produce SO’s with intermediate kinematic properties between fast and slow rotators 


that are difficult to explain with major mergers (Tapia et al. 20141. In conclusion the bulges of SO’s can be 


formed by means of a major merger of early disks or by a sequence of minor merger events, followed by a 
later disk rebuilding of the left-over gas and stripped stars (Somerville & Primack [[T999l. 


The data presented here are consistent with this scenario, although there are only hints that this was the 
right story for the SO in clusters. Direct evidences exist that ram pressure stripping and thermal evaporation 
could rapidly quench the SE in Spiral galaxies, while internal secular evolution can re-distribute the matter 
within them. 
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